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Abstract 


The Bose gas in an external potential is studied by means of the local density 
approximation. An analytical result is derived for the dependence of the critical 
temperature of Bose-Einstein condensation on the mutual interaction in a generic 
power-law potential. 
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In an ideal Bose-Einstein gas, the zero momentum state can become macroscopically occupied, 
and the system then undergoes a phase transition — the Bose-Einstein condensation (BEC) fi}. 
The A-transition of liquid helium-II is often interpreted essentially as a result of BEC, the strongest 
degeneracy effect of a Boson system f]. For many years, it was considered hopeless to experimentally 
observe BEC in an atomic gas with weak interactions. With the development of techniques to 
trap and cool atoms, BEC was recently observed directly in dilute atomic vapors, M, p- The new 
experimental achievements have stimulated great interest in the theoretical study of inhomogeneous 
Bose gases. 

For BEC in the atom-trapping experiments, the critical temperature Tę can be altered by the 
spatially varying potential of trapping. Furthermore, the noninteracting Bose gas can lead to 
unphysical results. The interaction between atoms may have a significant effect on the critical 
temperature. There have been several investigations analyzing the dependence of the critical tem- 
perature on the trapping potential and weak interaction in the Bose gas|6, [m]. Here we shall derive 
an analytical result for the critical temperature of a trapped interacting Bose gas by means of the 
local density approximation. 

When the energy level spacing of the trapping potential is much smaller than kT = 87t , the 
local density approximation B, oi is a good approximation. Space may then be divided into 
small cells, and in each cell we may consider the trapping potential energy V(r) to be a constant. 
In this approximation for an ideal Bose gas, the local density p(r) of particles in the gaseous phase 
(i.e., without BEC) is 

plr) =A? gsja(€), (1) 
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where 


€ = zexp[—GV(r)], z= PH, lz) = Site, (3) 


and y is the chemical potential. Equivalently, the density of states at energy € isf] 
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We may now extend the argument of BEC for rigid walls[i] to the case of the spatially varying 
trapping potential V(r). From Eq. M, the total number of particles can be written as 


N=N+ | për, (5) 


where No is the number of particles in the ground state. Let us consider a generic power-law 
potential as in Ref. ral 
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Since g3/2(€) is a bounded, positive, monotonically increasing function of €, for the potential (a), 
the integral in Eq. () reaches its maximum value at fugacity z = 1 or u = 0 . From Eqs. (8) and 


the critical temperature Te, at which the ground state begins to take on macroscopic values, is 
determined by 


V(r) =€ + €2 + €3 
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The definition of the Gamma function 
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For the specific choice a) of the potential, expression (a) then becomes 
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where 
I(p,1,q) = (pla) PL/p)P(A/D0 1/4), (11) 
and ¢(v 14” is the Riemann zeta function and 7 = 1/p + 1/l + 1/q + 1/2 . Noticing the 


ae a of the thermal wavelength on the temperature, we finally obtain 
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From Eqs. and ([0) we can obtain the fraction of condensation at a temperature T below Te 
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where No is the number of Bose particles in the BEC state. 
For the harmonic potential 
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the critical temperature and fraction of condensation reduce to 


where we have made use of 


¢(3) =1.202,  I(2,2,2) = n°/?/8. 
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The results (15) and (fE) were first obtained in Ref. and have been verified in a recent 


experiment 


(12). 
a for the generic power-law potential may also be derived in terms of Eq. H). In 


Ref. it is given that 


h3 N e,/P ef el? 


1/(n+1) 
2n(2m)3/? Lı LoL3 F(p,L, | 


kT, = | 


where 


co gn 
awn) = f aoa” 


1 i , 
F(p,l,q) =8 p dx (1 — eSa f dz (1 — ceria] ff dx (1— 2) '/2] . 
0 0 0 


In fact, we may integrate both Q and F to obtain 
Q(n) =f area = Cn +1) P+ 1), 
j=l 


F(p,l,q) = 4v7 I(p,1,q)/T(n + 1), 


where in deriving the last equation we have used 
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It can then be verified that Eq. (LẸ) coincides with (LJ) as it should. 
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Along similar lines we may analyze the critical temperature of a trapped interacting Bose gas. 
In Ref. it has been given that when the s-wave scattering length |a| < A, in the local density 
approximation the local density of the gaseous phase as given in Eq. m) is replaced by 


plr) = A393 /9(E), (24) 


with 7 
€ = zexp[-GV(r) — 4ad’p(r)], (25) 


where the term containing the scattering length a describes the mutual interaction of atoms. 


From Eqs. (24) and (25) we have 
p(0) = 8 93/2[z exp(—4ar° p(0)]. (26) 
Thus, at the critical temperature the fugacity is 
Ze = exp[4aX2p¢(0)] (27) 


where - 
pe(0) = AS °g3/2(1): (28) 
For this ze Eq. (RA) becomes 
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where we have used d 
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and kept only the lowest order in a/A. The determination of critical temperature T, for the 
interacting Bose gas involves integration of p(r) in 9). After some calculations similar to those 
in the derivation of (1d) from (A), we find 
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We finally obtain the critical temperature T, 
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Comparing Eq. B3) with Eq. (E) for the critical temperature Te of the trapped ideal Bose gas, we 
can derive the relation between the critical temperatures of BEC for ideal and nonideal Bose gases 
in the power-law potential 
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From this analytical relation between the critical temperature and the scattering length a, we see 
that the critical temperature T, decreases for a > 0, and increases for a < 0, compared with T, of 
the ideal Bose gas. This agrees with Ref. a. 

We may also obtain an expression similar to the case of an ideal Bose gas for the condensate 
fraction at a temperature T < T, to be 


Mi eu +4 esi xt M12) ow i - (=) "] (36) 


Te Te 


For the harmonic potential (4), Eq. (BG) reduces to 
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which has been obtained in Ref. [i]. 

The main effect of an external potential is to concentrate Bose particles. A repulsive mutual in- 
teraction between particles will weaken the effect of the external potential in concentrating particles 
around the center of the potential, hence lower the critical temperature. 
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